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Abstract 

Recent investigations by Bender and Boettcher (Phys. Rev. Lett 80, 5243 
(1998)) and Mezincescu ( J. Phys. A. 33, 4911 (2000)) have argued that the 
discrete spectrum of the non-hermitian potential V{x) = —ix'^ should be real. 
We give further evidence for this through a novel formulation which transforms 
the general one dimensional Schrodinger equation (with complex potential) 
into a fourth order linear differential equation for |^'(x)p. This permits the 
application of the Eigenvalue Moment Method, developed by Handy, Bessis, 
and coworkers (Phys. Rev. Lett. 55, 931 (1985);60, 253 (1988a,b)), yielding 
rapidly converging lower and upper bounds to the low lying discrete state 
energies. We adapt this formalism to the pure imaginary cubic potential, 
generating tight bounds for the first five discrete state energy levels. 
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I. INTRODUCTION 



In the recent work by Bender and Boettcher (1998) they conjectured that certain VT 
invariant systems should have real discrete spectra. Various examples were presented, in- 
cluding the —ix^ potential. The interest in such systems has increased, particulary through 
the more recent work of Bender et al (1999), Bender et al (2000), Bender and Wang (2001), 
Caliceti (2000), Delabaere and Pham (1999), Delabaere and Trinh (2000), Levai and Znojil 
(2000), Mezincescu (2000, 2001), Shin (2000), and Znojil (2000). 

We present a radically new way of attacking such problems. Although the results pre- 
sented here combine rigorous mathematical theorems and their numerical implementation, 
it should also be possible to develop them purely within an algebraic context, and confirm 
that the —ix^ potential can only have real discrete spectra. This particular approach is un- 
der investigation, and the results will be presented elsewhere. However, we have been able 
to implement the procedure discussed below, numerically, for the case of complex energies, 
E, and find no evidence for such discrete states (for moderate energy values). The details 
of this will be communicated in a forthcoming work focusing on the ix^ -\- iax potential 
studied by Delabaere and Trinh (2000). Our principal objective in this communication is to 
emphasize the importance of positivity as a quantization condition, within the appropriate 
(moment based) representation. 

Our starting point is the observation that the one dimensional Schrodinger equation (on 
the real hne), 

- dl^{x) + V{x)^{x) = E^{x), (1) 

for complex potentials, V — Vr + iVj, and real energies, ImE — 0, can be transformed into 
a fourth order, hnear differential equation for S{x) = |^'(x)p : 

+ {m + 2{^))s^0, (2) 



2 



where S^^^ = d^^S. This equation assumes that the eigenenergy, E, is real. We derive it in 
the next section. 

We could also assume that E is complex and incorporate its imaginary part into Vj. 
Since our objective is to show, numerically, that the conjecture that E is real is a viable one, 
we restrict our considerations to this case only, here. The method presented in this work is 
so powerful (both theoretically and numerically) that if the discrete state is not purely real, 
then it will be detected, at some sufficiently high calculation order. 

The above fourth order differential equation can be generalized to include any complex 
contour in the complex plane. However, for the particular problem considered here, we have 
only focused on the simplest representation for S{x), as given by Eq.(2). 

If the potential is real, Vi — 0, then Handy et al (1987a,b; 1988c) have shown that S{x) 
satisfies a third order differential equation. This is easy to see from the above by simply 
taking Vj — > 0, and recognizing that Eq.(2) becomes the total derivative of the third order 
equation 

- Is^^'Hx) + 2V{x)S^^\x) + V'{x)S{x) = 2ES^^\x). (3) 

The importance of converting the discrete state problem into the nonnegative S{x) repre- 
sentation is that for rational fraction complex potentials, one can then exploit the Eigenvalue 
Moment Method (EMM) of Handy, Bessis, and coworkers (1985,1988a,b), enabling the gen- 
eration of converging lower and upper bounds for the low lying discrete states. 

For rational fraction potentials, Eq.(2) can be transformed into a moment equation 
involving the Hamburger moments 

/+0O 
dx x^S{x), (4) 
-oo 

p > 0. The Moment Equation (ME) takes on the form 

ma 

lJip^Y.^pAE)y^i. (5) 

p > 0, where the energy dependent coefficients arc easily obtained, and satisfy (i.e. "ini- 
tialization conditions"): M^^^^j — ^luiii ^ ^1,2 < '^s- The missing moments, 

3 



{a*€|0 < ^ < rUs}, are to be considered as independent variables. The missing moment 
order, m^, is problem dependent. 

The homogeneous nature of the Schrodinger equation requires the imposition of an ap- 
propriate normalization condition. Although this requires some care, usually, a convenient 
choice is to take 

rUs 

E/^P = 1- (6) 
Solving for //q, and substituting into the ME relation, gives 

lJip = Y.M^,,{E)jle, (7) 

where 

ri, ^ = 

i^i = (8) 

and 

MpAE)^\ , , , , • (9) 

Prom the Hankel-Hadamard (HH) positivity theorems (Shohat and Tamarkin (1963)), 
the Hamburger moments must satisfy the conditions dx {Yj'j=QCjX^YS{x) > 0, for all 
C's and J > 0. These become the quadratic form expressions 

J 

E Cj,iij,+j,Cj, > 0. (10) 

ii,2=o 

In terms of the (unconstrained) normalized /x's this becomes 

nis J 

y: Cj,Mj,+j,ae)Cj,) > 0, (11) 

e=o ii,2=o 

which defines the linear programming equations (Chvatal (1983)): 

J2Ai{C,E)f,,<B{C,E), (12) 

e=i 

for all possible C's (except those identically zero), where 



MC,E) = -( E C,,M,,+,,,(i?)C,,), (13) 
h,2=o 

and 

BiC,E) = ( £ C,,M,,+,„o(i?)C,,). (14) 

jl,2=0 

If at a given order, J, and arbitrary energy value, E, there exists a solution set to all 
of the above inequalities, then it must be convex. Through a linear programming 

based cutting procedure (Handy et al (1988a,b)), one can find optimal C's which (in a finite 
number of steps) establish the existence or nonexistence of U^^\ The energy values for which 
missing moment solution sets exist, define energy intervals, 

NiJ) 

Ee U[^S,4;i], if^'Vo, (15) 

n=0 

which become smaller as J increases, converging to the corresponding discrete state energy 
(which must always lie within the respective interval): 

45 < 4f < ■ ■ • < E,Hys^cal■,n < • • • < < ^i- (16) 

Through the EMM approach, we can easily generate the converging lower and upper bounds 
to the desired discrete state energy. 

We note that although the traditional Moment Problem theorems are concerned with 
uniqueness questions (i.e. is there a unique function with the moments fXp satisyfing the HH 
positivity conditions ?), within the context of physical systems such issues are usually in- 
consequential. This is because the very nature of the ME relation will guarantee uniqueness. 
That is, our moments are associated with an underyling differential equation with unique 
physical solutions. 
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II. DERIVING THE POSITIVITY EQUATION FOR S{X) 

We derive Eq.(2) as follows. First, multiply the Schrodinger equation {E real) 

- ^*{x)^"{x) + V{x)S{x) = ES{x). 
The complex conjugate becomes 

- ^{x)^*"{x) + V*{x)S{x) = ES{x). 
Adding both expressions, and using ^*^" = (^*^'y _ |^'|2^ yields 

- [S" - 2|*f ] + 2VrS = 2ES. 
This in turn becomes (upon differentiating) 

- S'" + 2(|*f )' + 2{VrS)' = 2ES'. 
If we subtract Eq.(18) from Eq.(17), then 

=2iViS. 

Returning to the Schrodinger equation, we multiply both sides by ^ 
The complex conjugate is 

Substituting V = Vji + iVj, we add both expressions (and divide by iVj): 

_ J- + ^ + f^^*' _ = E— 

iVi iVi ^ ^ iVi 

Differentiating with respect to x, and substituting Eq.(21) yields 

^ iVi ^ ^ iVi ^ ^iVi^ 
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Upon dividing Eq.(19) by iVj, and differentiating, we obtain 




Finally, we substitute Eq.(25) for the second term in Eq.(26), obtaining a fourth order linear 
differential equation for S: 

or 

which becomes Eq.(2). 

The positivity differential representation in Eq.(2) is a fourth order linear differential 
equation, with four independent solutions, for any E. Within the EMM formalism, it is 
important to prove that the physical solution is the only one which is both nonnegative 
{S{x) > 0) and bounded, with finite moments (i.e. S{x) is in L^). We can prove this for 
Eq.(2). 

For any real energy variable value, E e let ^i(a;) and ^2(2^) denote the two indepen- 
dent solutions to the Schrodinger equation. The expression S{x) = \a'i/i{x) + /3^2(2;)P = 
lap X |*i(a:)p + |/3|2 x \-^2{x)\'^ + ap*'^fi{x)'^^{x) + a*p'^l{x)-^2{x), then becomes a solution 
to Eq.(2). So too are 1^^1(0;) P and 1^2(3:^) P- Accordingly, since a and (3 are arbitrary, and 
^i(x) and ^2(2^) are complex, the configurations ^'i(a;)\l'2(a;) and '^l{x)'^2{x) are indepen- 
dent (complex) solutions to Eq.(2) as well. 

From low order JWKB asymptotic analysis (Bender and Orszag (1978)), in either asymp- 
totic direction (x — > ±cxd), one of the semiclassical modes will be exponentially increasing, 
while the other is exponentially decreasing. Therefore it becomes clear that the only possible 
nonnegative and bounded S{x) configuration is that corresponding to the physical solutions. 
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III. THE -IX^ POTENTIAL 
The positivity differential equation for the V{x) — —ix^ potential is (i.e. Vr — 0,Vj — 

x-^S^^^ (x) - 3x-^S^^^ (x) + AEx-^S^^^ {x) - UEx'^S^^^ (x) - Ax^S{x) = 0. (29) 
Multiplying both sides by x^^^, and integrating over 3?, produces the ME relation 

4//j,+7 = {p + A)p{p - 1) (p - 2)/ip_3 + AEp{p + 4)^p_i, (30) 

for p>0. 

The moment equation separates into two relations, one for the odd moments, the other 
for the even moments. Assuming that the discrete states are nondegenerate and have real 
eigenenergies, we have: 

^f*{-x) = ^(x), (31) 

and 

S{-x) = ^*{-x)^{-x) = ^{x)^*{x) = S{x). (32) 

Thus, the physical 5'(a;)'s are symmetric, and the odd order moments are zero. 
The even order Hamburger moments 

At2p = Up, (33) 

correspond to the Stieltjes moments, 

POO 

up= / dyy'^Tiy), (34) 

J 



of the function 



-i(v) - (35) 



The corresponding Stieltjes moment equation for the —ix^ potential becomes (i.e. substitute 
p = 2p+l in Eq.(30)) 

Au,+4 = (2p + 5)(2p + l)(2p)(2p - l)up_, + AE{2p + l)(2p + 5)up, (36) 

for p > 0. This is an = 3 order problem. One can convert this into the form in Eq.(5) 
(i.e. Up = Yl^=o Mp^£{E)ui), where the M coefficients satisfy Eq.(36), with respect to the 
first index {p), as well as the initial conditions previously identified. 

One convenient feature about the Stieltjes representation is that the normalization con- 
dition 

E^^ = l' (37) 

e=o 

involves nonnegative moments. 

From the Stieltjes moment problem (Shohat and Tamarkin (1963)) we know that the 
counterpart to Eq.(lO) is 

h,2=o 

for cr = 0, 1. Accordingly, the necessary linear programming equations to consider are 

rris 

J2MC,E;a)<B{C\E;a), (39) 
e=i 

where 

MC,E;a) = -( E C,,M,+,,+,„,(E)C7,,), (40) 

jl,2=0 

and 

B{C,E;a) = ( E Cj,Ma+n+nAE)Cj,). (41) 

ii,2=o 

The numerical implementation of the EMM procedure yields the excellent results quoted 
in Tables I - V. Our results are in agreement with those of Bender and Boettcher (1998), as 
well as those of Handy, Khan, and Wang (2000). We indicate the maximum moment order 
generated, Pmax, through the ME relation. 
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Since our results are based on equations that explicitly assume E is real, and the EMM 
procedure is very stable and highly accurate (as evidenced through the tightness of its 
bounds), any imaginary part to the discrete state energy would reveal itself through some 
anomalous behavior in the generated bounds. That is, at some order Pmax, no feasible 
energy interval would survive (i.e. U^^ — 0, for all E). This is never observed, to the order 
indicated. As such, our analysis strongly supports the reality of the (low lying) discrete 
state spectrum for the —ix^ potential. 



10 



TABLES 



TABLE I. Bounds for the Ground State Energy of the —ix^ Potential 



max El-0 Eu-o 

10 .825 1.405 

20 1.15619 1.15645 

30 1.1562669 1.1562672 

40 1.1562670718 1.1562670721 

50 1.156267071988016 1.156267071988161 

60 1.15626707198811324 1.15626707198811335 



TABLE IL Bounds for the First Excited State Energy of the —ix^ Potential 



maa; -Ez,;! Eu-1 

20 4.1056 4.1168 

30 4.109225 4.109236 

40 4.1092287509 4.1092287578 

50 4.109228752806 4.109228752812 
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TABLE III. Bounds for the Second Excited State Energy of the —ix^ Potential 



moa; -E'l;2 -E'c/;2 

20 7.420 7.594 

30 7.56213 7.56242 

40 7.562273794 7.562273999 

50 7.5622738549 7.5622738551 



TABLE IV. Bounds for the Third Excited State Energy of the -ix^ Potential 

Pmax El-3 Eu-3 

30 11.3115 11.3159 

40 11.314418 11.314425 

50 11.314421818 11.314421824 



TABLE V. Bounds for the Fourth Excited State Energy of the —ix^ Potential 



30 15.20 15.80 

40 15.29145 15.29160 

50 15.29155366 15.29155380 

60 15.29155375037 15.29155375041 
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